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Abstract. Recent results on heavy flavor production at HERA are presented. New 
measurements of charm production are available in the region of Q 2 < 1000 GeV 2 
for D* and in the photoproduction limit, Q 2 — > 0, for D s . The results are compared 
to theoretical calculations and the charm proton structure function F% c (x,Q ) is ex- 
tracted. A new independent measurement of the bb cross section using a ^-lifetime tag 
is reported. 



I INTRODUCTION 

At the storage ring HERA positrons of 27.5 GeV collide with protons of 920 GeV 
resulting in a center of mass energy y/s = 320 GeV. The event kinematics of ep 
scattering is described by s = k + P, Q 2 = —(k — k') 2 , y = q ■ P/k ■ P, x = Q 2 /sy 
and W 2 p — (q + P) 2 = sy — Q 2 where k, k', P and q are the 4-vectors of the 
incoming/outgoing electron, incoming proton and exchanged photon, respectively. 



A Heavy Flavor Production Mechanism 



Heavy quark production predominantly occurs via the process of boson gluon 
fusion (BGF) sketched in fig. la. It is thus directly sensitive to the gluon distri- 
bution in the proton. At HERA the measurement of heavy flavor processes is a 
powerful means to study perturbative QCD down to values of x > 10 -4 . In par- 
ticular it can be tested if the parton distributions are universal, i.e. independent 
of the nature of the hard scattering process. The presence of two large scales, the 
squared 4-momentum transfer Q 2 and the heavy quark masses m 2 b makes heavy 
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FIGURE 1. a) Charm Production via Bo son- Gluon- Fusion; b) Mass difference distribution 
AM = Mkttits — Mktt &s obtained from 110 pb^ 1 by ZEUS. 



flavor production a useful test bed for the study of perturbative QCD resummation 
techniques. 

Several schemes to calculate heavy quark cross sections in next-to-leading order 
(NLO) are available. The program HVQDIS [1] implements a fixed flavor scheme 
at fixed order a 2 , assuming that the proton is made up of three active quark flavors 
(uds). The parton distributions are obtained by the DGLAP equation [2]. The 
heavy quark pairs are produced at the perturbative level via BGF. The results 
of HVQDIS are expected to be most accurate for Q 2 ~ m 2 c . At large Q 2 terms 
of orders higher than a 2 contain \og(Q 2 /m 2 ) factors that can become large. The 
CCFM equation [3] approximates the parton cascade as emitting gluons in an an- 
gular ordered way. The cross section is calculated according to the ^-factorization 
theorem by convoluting an unintegrated gluon density with an off-shell BGF ma- 
trix element. A new hadron level Monte Carlo generator CASCADE [4] is available 
which implements the CCFM equation. 



ZEUS 1995-2000 Preliminary ZEUS 1995-2000 Preliminary 




log 10 (Q 2 ) (GeV 2 ) log 10 (x) 



FIGURE 2. Differential D* production cross sections as a function ofQ 2 (a) andx (b) compared 
to the NLO QCD calculation of HVQDIS with different parton density distributions. 
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FIGURE 3. a) Distribution of M(K + K^n ± ) for events inside the <f) mass range 
1.0115 < M{K + K~) < 1.0275 GeV as measured by ZEUS; b) M(K+K~) distribution for events 
inside the D s mass range. 



B Experimental Procedure 

Charm events are tagged by reconstructing _D*-mesons (HI [5], ZEUS [6,7]), 
D s mesons (ZEUS [8]) or semi-leptonic decays into electrons (ZEUS [9]). Visible 
cross sections of D-meson production are measured in the central region of detector 
acceptance, at > 1.5 GeV/c and \rj\ < 1.5, where pj_ is the transverse momentum 
of the D-meson and rj denotes the pseudorapidity. Measurements of the 66-cross 
section have been performed through semi-leptonic decays in muons (HI [10], ZEUS 
[11]) or electrons (ZEUS [12]) and, more recently, by exploiting the long lifetime of 
5-hadrons (HI [13]). 

Since the beginning of HERA operation in 1992 the luminosity has continuously 
increased, yielding a total integrated luminosity collected by each of the two collider 
experiments HI and ZEUS of more than 100 pb _1 . Fig. lb shows the D* — D° 
mass difference distribution measured by the ZEUS experiment. There are roughly 
27,000 events in the signal illustrating the excellent capabilities of HERA for the 
study of the charm production mechanism and also for charm meson spectroscopy 
studies [14]. In winter 2000/2001 the HERA accelerator is being upgraded. In the 
future it is expected to deliver a yearly integrated luminosity of 150 pb _1 . 



II OPEN CHARM PRODUCTION 

A D* production at high Q 2 

Fig. 2 shows the results of 82.6 pb -1 of D* data collected by the ZEUS experiment 
[6] in the region 10 < Q 2 < 1000 GeV 2 . The shaded band represents the theoretical 
calculation as implemented in the HVQDIS program. Here the charm quark mass 
is varied between 1.3 and 1.6 GeV. The input gluon distribution is derived from 
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FIGURE 4. Differential cross section (a) da/dp^ and (b) dajdr\ for D* (triangles) and D, 

(bullets) in comparison with NLO calculations (see text). 



results of different QCD fits to the inclusive proton structure function [15-17]. The 
Peterson fragmentation function [18] is used to describe the probability that a c- 
quark hadronizes into a D* meson. The theoretical predictions agree well with the 
data indicating that, within the statistics presently available, the BGF calculation 
in the fixed flavor scheme of order a 2 fully accounts for the charm cross section up 
to very high values of x and Q 2 . 



B D* and D s production in photoproduction 

Visible cross sections of D* mesons in photoproduction have been presented be- 
fore [19]. More recently a measurement of _D s -mesons was performed using the decay 
D s — > 07r and — > KK. Fig. 3 shows the signals of the reconstructed D s and the 
corresponding <fi [8]. The visible cross section is extracted in the region 12 > p®* > 3 
GeV and \r] D "\ < 1.5 yielding a%° x = (3.79 ± 0. 59 (stat)± 2 4 l(syst) ± OM(BR)) nb 
where the last error accounts for the uncertainty of the branching ratio. Fig. 4 
compares the measured visible production cross sections to NLO predictions [20]. 
The bold curves refer to a charm mass value m c = 1.5 GeV and the renormal- 
ization scale m# = a^J m 2 + p\ with a = 1. From fig. 4a it is evident that the 
normalization of the theoretical prediction is significantly too low even for ex- 
treme values m c = 1.2 GeV and a = 0.5 (thin lines). Furthermore the shape 
(fig.4b) is inconsistent with the data which are much higher in the forward di- 
rection 77 > 0. This behavior of the data could be due to an interaction be- 
tween the color charges of the c-quark and the proton remnant (beam drag) [21] 
which is not accounted for in the model. With the PYTHIA fragmentation the 
agreement is slightly better [22]. In the same kinematic region a result of a™ = 
(9.17 ±0.35(stat)1loig(syst)) nb is obtained for D* production. The ratio a^Ja'gi 
is determined as 0.41±0.07(s?/st)loQ5(stat) ±0.10(Si?). This result is compared to 
the ratio of fragmentation probabilities f(c -> D s )/f(c — >• D*) = 0.43 ± 0.04 ± 0.11 




FIGURE 5. Differential cross sections for D* production from the Kixix final state (bullets). 
Fig. d and e also show the cross section for the KA-k decay channel. 



as obtained from e + e~ experiments [23]. The strangeness suppression factor 7 S 
is the ratio of probabilities to create s or (u, d) quarks during the fragmentation 
process. In simulations based on the Lund String Model [24] 7 S can be adjusted. 
For HERA an optimal value of 0.27 ± 0.05 is found while for e + e~ experiments 
0.26 ± 0.03 had been obtained [23]. While direct measurements of the fragmen- 
tation function at HERA are not yet available, the agreement shows a consistent 
picture of charm fragmentation which tends to support the universality of charm 
fragmentation. 



C D* cross sections at small Q 2 

Visible D* cross sections at small values of Q 2 have been measured by both HI [5] 
and ZEUS [7] yielding <t™™ d . x = 8.31 ± 0.31(stat)t° i(syst) nb and a^ eD * x = 
8.37 ± 0Al(stat)±lil(syst)±° 6 3 g(theo) nb for slightly different kinematic regions 
(ZEUS: 1< Q 2 < 600 GeV 2 , y < 0.7, 1.5 < p ± < 15 GeV,M < 1.5; HI: 1 < Q 2 < 
100 GeV 2 , 0.05 < y < 0.7, p ± > 1.5 GeV,\r]\ < 1.5). The ZEUS data are shown 
in fig. 5 in comparison with HVQDIS. The open band corresponds to the standard 
Peterson fragmentation function. To account for potential beam-drag effects the 
Peterson fragmentation is replaced by that of the RAPGAP simulation [25] (shaded 
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FIGURE 6. Single differential inclusive cross section a(ep — > eD*X) as a function of a) Q 2 , 
b) x, c) r/. The expectations from HVQDIS (CASCADE) are indicated by the dark (light) shaded 
bands. 

band) leading to a somewhat improved description of the r] D * distribution. Fig. 6 
shows the HI data in comparison with the result from HVQDIS (dark shaded) and 
from CASCADE (light shaded). CASCADE describes the data fairly well in both 
normalization and shape while HVQDIS is lower and fails to describe the shape 
of the t)d* -distribution in particular in the forward direction. At small values of x 
CASCADE also provides a better description. 

Ill EXTRACTION OF Ff° 

The charm contribution to the proton structure function, F£ c , is related to the 
charm cross section in x and Q 2 as = + (1 ~ vY) ' Ff( x ,Q 2 )- The 

contribution from Fl is expected to be small and is neglected here. _F 2 CC is extracted 
from the measured visible cross section in each bin of x and Q 2 by extrapolation 
into the full phase space in rj and p± using the HVQDIS program. As shown in 
fig. 7 there is good agreement between HI and ZEUS and, in general, with the 
prediction obtained from a NLO QCD fit to the inclusive HI data. At the lowest 
x and Q 2 , a tendency towards a departure from the model is visible. Alternatively 
the CCFM model as implemented in CASCADE can be used for the extrapolation. 
Fig. 7b) shows the comparison between the HI data as extracted using CCFM and 
the CCFM expectation. There is good agreement also at the lowest values of x 
and Q 2 . A comparison of F 2 C0 with the inclusive proton structure function F 2 shows 
that the contribution from charm to the inclusive cross section is larger than 10% 
everywhere in the region of Q 2 > 1 GeV 2 and x < 10 -3 and can be as large as 25% 
at large values of Q 2 . 

IV OPEN BEAUTY PRODUCTION 

In contrast to charm quarks, 6-quarks are much heavier, thus providing a harder 
scale for perturbative calculations. At the same time the total bb cross section is ex- 



a) 



in the NLO DGLAP scheme 



HI NLO QCD 

fit to F, 



0.50 



0.25 



• HI Prelim. □ ZEUS 

V ZEUS e" prelim. 



0.00 



O CI 



0.40 



0.00 



0.40 



0.00 



Q = 1.5 GeV" 
(1.8) 
(2.5) 



_3e 



f-R — i — h 




60 GcV" 



3 .5 GeV" 
(4) 



18 GeV" 

(20) 



H 1 1 1- 



130 GeT 
\ , (100) 




w 

O 

< 
U 

< 
U 



b) 

0.40 
0.20 
0.00 
0.25 
0.00 

i 

0.40 
0.00 



F 2 c in the CCFM scheme 



A HI Preliminary 



CCFM 



1 1 1 1 

Q 2 =1.5GeV 2 ■ 


i 


i i i 

3.5 GeV 2 ■ 




i 




: — 1 — 1 — 1 — 1 — : 

6.5 GeV 2 ■ 


; h- 


— 1 1 1 ; 

12 GeV 2 ■ 


: \ 

i 

1 1 1 1 1 '■ 


'■ h- 


— 1 1 1 '■ 


25 GeV 2 




60 GeV 2 . 


■ N ; 

i i i i 


i 


k : 
i i i 



-5 -4 -3 -2 -1 -5 -4 -3 -2 -1 -5 -4 -3 -2 -1 

logx 



-1 



-5 -4 -3 -2 -1 -5 -4 -3 -2 

logx 

FIGURE 7. T/ie charm contribution to the proton structure function F£ c as derived a) using the 
HVQDIS program (HI, ZEUS) b) using the CCFM model as implemented in CASCADE (HI). 
Also shown are the predictions for F^ varying the charm quark mass between 1.3 and 1.5 GeV. 



pected to be smaller than that of charm by a factor of 200. The first measurements 
on the beauty production cross section in ep-scattering at HERA were performed 
using semi-leptonic decays [10,12,11]. In these analyses events with two or more 
jets are selected if they have at least one identified muon (electron) with a p± of 2 
GeV (1.6 GeV) in the central region of the HI (ZEUS) detector respectively. The 
lepton is assigned to one of the two jets of the event. To distinguish 6-events from 
those with charm or fake leptons and to extract the bb cross section the transverse 
momentum p™ 1 of the lepton relative to the jet is used (fig. 8a). Both measurements 
(fig.9) are above the NLO expectations. 

Recently a new independent measurement has been presented by HI [13]. The 
analysis uses the central silicon vertex detector to improve the track resolution. 
Events with 6-mesons are tagged by their long lifetime using the impact parameter 
between the muon track and the beam in the transverse plane. From the impact 
parameter distribution (fig. 8b) the 66-cross section is extracted. A log-likelihood 
fit yields a value of cr ep ^ eh i x = (176 ± 30(stat) ± 29(syst)) pb which confirms the 
previous results. To maximize the 6-contribution in the event sample and to reduce 
the systematics the two methods are combined in a likelihood fit (fig. 10), yielding a 
cross section of 170±25 pb. Like the previous measurements this result is somewhat 
higher than the NLO QCD expectation [20] which gives 104 ± 17 pb. 




V SUMMARY 



Since 1992 HERA has been collecting a wealth of data at continuously increasing 
luminosity. New results with charm and beauty data have been presented here. In 
the charm sector precision tests of QCD are being performed by confronting the 
data in a variety of channels (D*, D s and semi-leptonic decays) with predictions 
from NLO calculations. The process of boson gluon fusion has been shown to be by 
far the dominant contribution to the production of charm quarks in the wide kine- 
matic range between the photoproduction regime, Q 2 — > 0, and Q 2 < 1000 GeV 2 . 
However, in photoproduction, presently, the models do not fully reproduce the 
data in normalization and shape. In the region of low Q 2 (1 < Q 2 < 600GeV 2 ) the 
charm contribution to the proton structure function F 2 CC has been extracted. The 
data have been compared to calculations using different parton evolution schemes 
and there are indications that the CCFM equation gives a better description of 
the data in the region of the smallest values of x. In the beauty sector a new sys- 
tematically independent measurement using the lifetime tag (impact parameter) 
of semi-muonic events has been presented. The result confirms the previous mea- 
surements from HI and ZEUS which are somewhat larger than expected by NLO 
calculations. 
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FIGURE 9. a) Extrapolated b photoproduction cross section as measured by the ZEUS Experi- 
ment for Pt > p™ m = 5 GeV, 0.2 < y < 0.8 and \r] b \ < 2 compared with NLO QCD predictions 
plotted as a function of p^ m . b) The total photoproduction cross section, a(^p — > bbX) measured 
by HI. The horizontal error bar represents the range of the measurement. The curve shows the 
expectation from the FMNR NLO calculation using different structure functions. 
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FIGURE 10. a) Impact parameter distribution for muon candidates with p r ^ 1 > 2GeV and 
p r ^ 1 distribution for muon candidates with S > 500/.im, with estimated contributions. 
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